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Abstract 


A united approach of the large-scale structure of a closed universe and the lo- 
cal spherically symmetric gravitational field is given by supposing an appropriate 
boundary condition. The general feature of the model obtained are the following. 
The universe is approximately homogeneous and isotropic on the average on large 
scale and is expanding at present, as described by the standard model; while locally, 
the small exterior region of a star started long ago to contract, as expected by the 


gravitational collapse theory. 
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1 Introduction 


In the traditional approaches of the universe model and the local gravitation problem, 
the two subjects are alone treated respectively. In various universe models people only 
consider the average structure on large scale; the distribution of matter is imagined as 
homogenized or smoothed. A galaxy is nothing but a perfect fluid element. On the other 
hand, when people investigate local gravitation problem such as the spherically symmet- 
rical gravitational field, they often regard the region as asymptotically flat spacetime, as 
if there is no remain matter in the universe. True, such a treatment has its reasons: (i) 
the model is simple, and it does not impair the knowledge of the large-scale structure of 
the universe to ignore its minor details. (ii) For a local gravitation problem, the Birkhoff’s 
theorem [I] seems to be the basis for leaving the remaining matter of the universe out 
of consideration. However, no matter how reasonable it is apparently, the asymptotically 
flat behavior of the solution of local gravitational field does not be in tune with a closed 
universe model, while the real universe most likely is closed. 

The cosmological redshift shows that the distances between galaxies are increasing, 
that is, the universe is expanding on the scale of the space between galaxies at present. 
On the other hand, the existing main form of the cosmic matter is the celestial body, some 
of which are very compact. This shows that in some very small areas of the expanding 
universe the space has begun to contract long ago. To our knowledge, there is yet no model 
to treat the two problems concurrently at present. In the present article we try to give, 
by means of the investigation a spherically symmetric gravitational field in an expanding 
closed universe, a united approach to the universe model and a local gravitational problem. 

In order to make it practical, we shall adopt a compromise proposal, which can be 
called “homogeneous by areas model” of matter distribution, that is, the universe can be 
divided into a few spherically symmetric parts with the same symmetric center, and in 


each part the matter density is a constant. 


2 The Gravitational Field Equations and Their Gen- 


eral Solutions 


Imagine in a closed Friedmann universe a local dense-matter area appears owing to a 


certain kind of mechanism. For simplicity,we assume that its density is a constant. Of 


course, around the dense-matter area a spherically symmetric rare-matter or vacuum area 
may appear. Such a matter distribution model is our start point. In this case, the 3- 
dimension spacelike hypersurface, t = constant, is a rotational hypersurface embedded in 


a 4-dimension flat spacetime. In comoving coordinates the line element can be written 
ds? = dt? — U (y, td? — V (Y, t) sin? W(d6? + sin? 0d¢”) (1) 


where0 << w<7,0<0<27,0<@< 2r. The relation between 4-dimension spherical 


coordinates R, Y, 0 and @ and the Decare coordinates x, y, z and w are 


x = Rsinysinð coso 
y = Rsinwsinésing 
z = Rsinwcosd 


w = Rcosy (2) 
One often introduces 
r = siny (3) 
and rewrites (I) as the form 
ds? = dt? — A(r,t)dr? — B(r, t) (d8? + sin? 0d¢”) (4) 


which is a well-known form [3]. 

However, it is necessary to emphasize that r is not a good coordinate owing to non- 
monodromy of the transformation (8). The same value of r corresponds to two values of 
p, i.e., two positions in the universe. This means that A(r,t) and B(r,t) should have two 
set physical solutions in principle, which corresponds to the two areas 0 < Y < 7/2 and 
m/2< < r, respectively. 


Consider the zero pressure perfect-fluid model 
Le = alle” (5) 
where U” = U?’ = U? = 0, Ut = 1. The Einstein’s equations are [I] 
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From the last equation of (6), it is easy to find that 
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here K = K(r) is an arbitrary function of r. Using the firt three equations of (6) in (@ 
gives that - 
BB? 4Kr’? 
4 = 8 
Set 
B = S(r, t)r’ (9) 
equation (8) becomes 
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from which we find se 
s? = oe (10) 


here F = F(r) is an arbitrary function of r. The solution of equation (IQ) is cycloid-like 


S F(1 — cosa) 
C+t = F(a-sina)/VK (11) 


where C = C(r) is an arbitrary function. Substituting equations (J), (Q) and (Ui) into 
(G) gives the expression of density 


3 (F Kr)! 


~ mG Sr (12) 
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3 Homogeneous by Areas Model of the Distribution 
of Matter 


Where we investigate the problem of spherically symmetric gravitational collapse of a 
celestial body in a closed universe, the ” homogeneous by areas model” above is a suitable 


model. This model finds expression in the form of the density, that is 


p(t) OSY <p (area I) 
pír, t)=4 0 pı <Y<y<Tr/2 (area IT) (13) 
Pult) Y< YIT (areal IT) 
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where Yı = arcsin r1, Ys = arcsin rə. The joint conditions on the boundaries Y% = yı and 
Wp = p2 is taken as that the metric g,,, are continuous on these boundaries. Substituting 


the expression into equation (12) and integrating gives 


Kr = ££ p, F?(1 — cos ar)’ (area T) 
Ky = A/(Fur®) (area IT) (14) 


Kor = 8 oF? (1 — COS arr)’ (areal IT) 


where A is the integral constant for the area JI; for the area J and the area JIT, the 
integral constants have been taken as zero so that both AK; and Krzr are finite at r = 0 


(i.e., Y = 0 and Y = r). The joint conditions are 


Br(ri,t) = Br(ri,#), Brr(r2,t) = Brr(r2,t) (15) 


The value of the parameter a in equation (QI) decides that the scale factor S(r, t) is 


expansive or contractive. As an additional condition, we demand that 
ar(rı, t) =ar(71,t),  Arr(r2,t) = arr1(T2, t) (16) 


The expressions (15) and (6) are determining solution conditions to F(r), K(r) and 
C(r). From (15) we can infer that S(r,,t) is continuous on the boundaries and then from 
the first of and (16) we can know 


Fy(r1) = Fri(mi), Frr(r2) = Frrr(r2) (17) 
The second of (QI) can be written as 
t—ty = F; (as — sin a;) — (ajo — sin aio)| [V Ki (18) 


here i = I,II, ITI and ajo = a;(r, to). Equations (6), (7) and (18) imply 


Kr(rı) = Kr(rı), Kr(r2) = Krrr(r2) (19) 
Substituting C6), C and C9) into M) gives 
4 3 
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From the first and the last of (I4 it is also easy to see that 


ps(1—cosayz)® = pso(1 — cos aro)? 


Pall — cosarr)? = pyo(1 — cos azz)? (20) 


In general, these conditions can yet not determine F; and K; (or ajo) completely 
because there still exists a certain freedom of coordinate transformation in comoving 
coordinates. In this sense, to select the function form of F; and ajg is just to select 
coordinates system. 

Now we construct a set of concrete solutions from the following considerations: (i) In 
the case Y2 = Yı (Ps = Pu, see (20)), the universe should be homogeneous and isotropic. 
(ii) In the case Y >> Y1, arro and &zrro Should be approximate constants and so are Fyr and 
Frrr. (iii) In the case Y & Yz both azo and &zro should be larger than 7, so that the scale 


factor in this area has been contracted, which reflects gravitational collapse of a celestial 


body. The first condition implies œo = Qy79 = Qro = Constant, and Fr = Fy; = Frrr 


constant. We try to choose the following forms of ajo and F; 


Qio = B+ ra~ "1 be t/y I2 (21) 
‘ T2 + Ty 

F;=a— eek be P/V bib2 (22) 
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where 3, 6, a and b are constants. Obviously, the last of (15) is satisfied because of the 
conditions (£6), (7) and (19). In order to make the first of (IS) held, we differentiate 
([8) with respect to r and substitute into ([4). So doing, we obtain 
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sin &zr) — (azo — sin arro)| 


(1 — cosazr)ah, = (1—cosaro)Qrry + 


3 sin arrr 
(1 — cos arrr)ahr = fa — cos ayzr70) + : 
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2(1 — COS Q]TIO 
sin arr) = (aærrro — sin arr0)| bairn (23) 


The continuity of œ on the boundaries gives 
(Fori 
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i i 
Ee „o l a arr : (24) 
or equivalent forms 
bwi 1 in ee dw, sin( 8 + dw) 1 
(a — bwi) vipa f1-r? "i [1 — cos(3 + dur)] Vides V1- r7 
bwo 1 i 1 _ dw sin( 8 + dw) 1 (25) 
(a — bw) Vide y1- r} r2 f — cos(8 + õwa)| Vida V1 -73 
where 
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If we assign any two of the four constants, say 8 and a (initial conditions), the other two 


can be determined by (25), and a set of uniquely definite solutions is thus found. 


4 Disussion 


From (27) and (22) it is easy to see that if rə = rı then aj = 8 and F; = a, which exactly 
is the homogeneous and isotropic Friedmann universe model. In general, where w > yı 
and Y >> we, we have ayy % 8 and Frrr © a, which shows that the spacetime far away 
from the star tends to become homogeneous and isotropic. In the place where Y < we 
We have 


T2 — rı 
N —<—— ó 26 
aro a ET (26) 
T2 — ri 
Fy rea- b 27 
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The existing observational data show that 8 ~ 7/2. If (ra —11)/(r2 + rı) > 7/2, then 
the celestial body is collapsing. 

In the traditional treatment of the problem of the spherically symmetric gravitational 
collapse, the exterior region of the star is considered infinite and empty Ø 3) A. The 
surface of the star is only boundary. A puzzle arose out of this treatment, which we 
arrate as follows. Suppose the surface of a star with constant density is described by the 
equation 

ety? 2 


oS (28) 


in a certain coordinate system. Under the transformations 


eon, y =y; Z = azb (29) 
The surface equation becomes 


gl? + y’? + z2 = a2 (30) 


in the new coordinate system. Then which symmetry has the gravitational field, axisym- 
metric or spherically symmetric? In our new treatment there is no such a puzzle owing 
to the existence of the second boundary. 

Birkhoff’s theorem states that the solutions of spherically symmetric and vacuum 
gravitational field equations must be Schwartzschild. Hence the metric of the area JJ 
should be able to become Schwartzschild form by means of a certain transformation. 


Such a transformation has been found, which is 


r = Fyr(l— cosarr) 
sa Seige eV ee) 
vV Frrr — GM sin arr — VGM(1+4+ cosarr) 
| Frrr — GM Frr— GM , 
aay ee aa Pnr pag nt 
GM F3 r3 l ' 
-j \/ reall aa (110 — sin or10)| dr (31) 
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